Search for: All records

Abstract A range of stellar explosions, including supernovae (SNe), tidal disruption events (TDE), and fast blue optical transients (FBOTs), can occur in dusty environments initially opaque to transients’ optical/UV light, becoming visible only once the dust is destroyed by transients’ rising luminosity. We present axisymmetric, time-dependent radiation transport simulations of dust-shrouded transients withAthena++and tabulated gray opacities, predicting the light curves of the dust-reprocessed infrared (IR) radiation. The luminosity and timescale of the IR light curve depend on whether the transient rises rapidly or slowly compared to the light-crossing time of the photosphere,t_lc. For slow-rising transients (t_rise ≫ t_lc) like SNe, the reprocessed IR radiation diffuses outward through the dust shell faster than the shell sublimates; the IR light curve therefore begins rising prior to the escape of UV/optical light, but peaks on a timescale ∼t_riseshorter than the transient duration. By contrast, for fast-rising transients (t_rise ≪ t_lc) such as FBOTs and some TDEs, the finite light-travel time results in the reprocessed radiation arriving as an “echo” lasting much longer than the transient itself. We explore the effects of the system geometry by considering a torus-shaped distribution of dust. The IR light curves seen by observers in the equatorial plane of the torus resemble those for a spherical dust shell, while polar observers see faster-rising, brighter, and shorter-lived emission. We successfully model the IR excess seen in AT2018cow as a dust echo, supporting the presence of an opaque dusty medium surrounding FBOTs prior to explosion. 

Abstract We expand the theoretical framework by O. Gottlieb et al., which connects binary merger populations with long and short binary gamma-ray bursts (lbGRBs and sbGRBs, respectively), incorporating kilonovae (KNe) as a key diagnostic tool. We show that lbGRBs, powered by massive accretion disks around black holes (BHs), should be accompanied by bright, red KNe. In contrast, sbGRBs—if also powered by BHs—would produce fainter, red KNe, potentially biasing against their detection. However, magnetized hypermassive neutron star (HMNS) remnants that precede BH formation can produce jets with power (P_NS ≈ 10⁵¹erg s⁻¹) and Lorentz factor (Γ > 10) likely compatible with sbGRB observations, and would result in distinctly bluer KNe, offering a pathway to identifying the sbGRB central engine. Recent modeling by J. C. Rastinejad et al. found luminous red KNe consistently accompany lbGRBs, supporting their origin in BH-massive disk systems, likely following a short-lived HMNS phase. The preferential association of sbGRBs with comparably luminous KNe argues against the BH engine hypothesis for sbGRBs, while the bluer hue of these KNe provides additional support for an HMNS-driven mechanism. Within this framework, BH–NS mergers likely contribute exclusively to the lbGRB population with red KNe. Our findings suggest that GW170817 may, in fact, have been an lbGRB to on-axis observers. Finally, we discuss major challenges faced by alternative lbGRB progenitor models, such as white dwarf–NS or white dwarf–BH mergers and accretion-induced collapse forming magnetars, which fail to align with observed GRB timescales, energies, and KN properties. 

Abstract Periodic collisions between a star on an inclined orbit around a supermassive black hole and its accretion disk offer a promising explanation for X-ray “quasiperiodic eruptions” (QPEs). Each passage through the disk midplane shocks and compresses gas ahead of the star, which subsequently re-expands above the disk as a quasi-spherical cloud. We present spherically symmetric Monte Carlo radiation transport simulations that follow the production of photons behind the radiation-mediated shock, Comptonization by hot electrons, and the eventual escape of the radiation through the expanding debris. Such 1D calculations are approximately justified for thin disks (scale-heighth ≲ few × R_⋆), through which the star of radiusR_⋆passes more quickly than the shocked gas can flow around the star. For collision speedsv_coll ≳ 0.15cand disk surface densities Σ ∼ 10³g cm⁻²characteristic of those encountered by stellar orbits consistent with QPE recurrence times, the predicted transient light curves exhibit peak luminosities ≳10⁴²erg s⁻¹and Comptonized quasi-thermal (Wien-like) spectra that peak at energieshν ∼ 100 eV, which is broadly consistent with QPE properties. For these conditions, gas and radiation are out of equilibrium, rendering the emission temperature harder than the blackbody value due to inefficient photon production behind the radiation-mediated shock. The predicted eruptions execute counterclockwise loops in hardness–luminosity space, qualitatively similar to QPE observations. Reproducing the observed eruption properties (duration, luminosity, temperature) requires a large radiusR_⋆ ≳ 10R_⊙, which may point to inflation of the star’s atmosphere from repeated collisions. 

Abstract We present nucleosynthesis and light-curve predictions for a new site of the rapid neutron capture process (r-process) from magnetar giant flares (GFs). Motivated by observations indicating baryon ejecta from GFs, J. Cehula et al. proposed that mass ejection occurs after a shock is driven into the magnetar crust during the GF. We confirm using nuclear reaction network calculations that these ejecta synthesize moderate yields of third-peakr-process nuclei and more substantial yields of lighterr-nuclei, while leaving a sizable abundance of free neutrons in the outermost fastest expanding ejecta layers. The finalr-process mass fraction and distribution are sensitive to the relative efficiencies ofα-capture andn-capture freeze-outs. We use our nucleosynthesis output in a semianalytic model to predict the light curves of novae breves, the transients following GFs powered by radioactive decay. For a baryonic ejecta mass similar to that inferred of the 2004 Galactic GF from SGR 1806-20, we predict a peak UV/optical luminosity of ∼10³⁹–10⁴⁰erg s⁻¹at ∼10–15 minutes, rendering such events potentially detectable to several Mpc following a gamma-ray trigger by wide-field transient monitors such as ULTRASAT/UVEX. The peak luminosity and timescale of the transient increase with the GF strength due to the larger ejecta mass. Although GFs likely contribute 1%–10% of the total Galacticr-process budget, their short delay-times relative to star formation make them an attractive source to enrich the earliest generations of stars. 

Abstract Collapsars—rapidly rotating stellar cores that form black holes—can power gamma-ray bursts and are proposed to be key contributors to the production of heavy elements in the Universe via the rapid neutron capture process (r-process). Previous neutrino-transport collapsar simulations have been unable to unbind neutron-rich material from the disk. However, these simulations have not included sufficiently strong magnetic fields and the black hole (BH), both of which are essential for launching mass outflows. We presentνh-amr, a novel neutrino-transport general relativistic magnetohydrodynamic (νGRMHD) code, which we use to perform the first 3D globalνGRMHD collapsar simulations. We find a self-consistent formation of a weakly magnetized dense accretion disk, which has sufficient time to neutronize. Eventually, substantial magnetic flux accumulates near the BH, becomes dynamically important, leads to a magnetically arrested disk (MAD), and unbinds some of the neutron-rich material. However, the strong flux also hinders accretion, lowers density, and increases neutrino-cooling timescale, which prevents further disk neutronization. Typical collapsar progenitors with mass accretion rates, $\dot{M}\sim 0.1-1M_{\odot }{\mathrm{s}}^{-1}$, do not produce significant neutron-rich (Y_e < 0.25) ejecta. However, we find that MADs at higher mass accretion rates, $\dot{M}\gtrsim \text{few}{M}_{\odot }{\mathrm{s}}^{-1}$(e.g., for more centrally concentrated progenitors), can unbindM_ej ≲ M_⊙of neutron-rich ejecta. The outflows inflate a shocked cocoon that mixes with the infalling neutron-poor stellar gas and raises the final outflowY_e; however, the finalr-process yield may be determined earlier at the point of neutron capture freeze-out. Future work will explore under what conditions more typical collapsar engines becomer-process factories. 

Abstract The origin of heavy elements synthesized through the rapid neutron capture process (r-process) has been an enduring mystery for over half a century. J. Cehula et al. recently showed that magnetar giant flares, among the brightest transients ever observed, can shock heat and eject neutron star crustal material at high velocity, achieving the requisite conditions for anr-process. A. Patel et al. confirmed anr-process in these ejecta using detailed nucleosynthesis calculations. Radioactive decay of the freshly synthesized nuclei releases a forest of gamma-ray lines, Doppler broadened by the high ejecta velocitiesv ≳ 0.1cinto a quasi-continuous spectrum peaking around 1 MeV. Here, we show that the predicted emission properties (light curve, fluence, and spectrum) match a previously unexplained hard gamma-ray signal seen in the aftermath of the famous 2004 December giant flare from the magnetar SGR 1806–20. This MeV emission component, rising to peak around 10 minutes after the initial spike before decaying away over the next few hours, is direct observational evidence for the synthesis of ∼10⁻⁶M_⊙ofr-process elements. The discovery of magnetar giant flares as confirmedr-process sites, contributing at least ∼1%–10% of the total Galactic abundances, has implications for the Galactic chemical evolution, especially at the earliest epochs probed by low-metallicity stars. It also implicates magnetars as potentially dominant sources of heavy cosmic rays. Characterization of ther-process emission from giant flares by resolving decay line features offers a compelling science case for NASA’s forthcoming COSI nuclear spectrometer, as well as next-generation MeV telescope missions. 

Abstract Although Cherenkov detectors of high-energy neutrinos in ice and water are often optimized to detect teraelectronvolt–petaelectronvolt neutrinos, they may also be sensitive to transient neutrino sources in the 1–100 GeV energy range. A wide variety of transient sources have been predicted to emit gigaelectronvolt neutrinos. In light of the upcoming IceCube Upgrade, which will extend the IceCube detector’s sensitivity down to a few gigaelectronvolts, as well as improve its angular resolution, we survey a variety of transient-source models and compare their predicted neutrino fluences to detector sensitivities, in particular those of IceCube-DeepCore and the IceCube Upgrade. We consider ranges of neutrino fluence from transients powered by nonrelativistic shocks, such as novae, supernovae, fast blue optical transients, and tidal disruption events. We also consider fast radio bursts and relativistic outflows of high- and low-luminosity gamma-ray bursts. Our study sheds light on the prospects of observing gigaelectronvolt transients with existing and upcoming neutrino facilities. 

Abstract Some hydrogen-rich core-collapse supernovae (SNeIIP) exhibit evidence of a sustained energy source powering their light curves, resulting in a brighter and/or longer-lasting hydrogen recombination plateau phase. We present a semi-analytic SNIIP light-curve model that accounts for the effects of an arbitrary internal heating source, considering as special cases⁵⁶Ni/⁵⁶Co decay, a central engine (magnetar or accreting compact object), and shock interaction with a dense circumstellar disk. While a sustained internal power source can boost the plateau luminosity commensurate with the magnitude of the power, the duration of the recombination plateau can typically be increased by at most a factor of ∼2–3 compared to the zero-heating case. For a given ejecta mass and initial kinetic energy, the longest plateau duration is achieved for a constant heating rate at the highest magnitude that does not appreciably accelerate the ejecta. This finding has implications for the minimum ejecta mass required to explain particularly long-lasting SNe, such as iPTF14hls, and for confidently identifying rare explosions of the most massive hydrogen-rich (e.g., Population III) stars. We present a number of analytic estimates that elucidate the key features of the detailed model. 

Abstract Superluminous supernovae (SLSNe) radiate ≳10–100 times more energy than ordinary stellar explosions, implicating a novel power source behind these enigmatic events. One frequently discussed source, particularly for hydrogen-poor (Type I) SLSNe, is a central engine such as a millisecond magnetar or accreting black hole. Both black hole and magnetar engines are expected to channel a fraction of their luminosity into a collimated relativistic jet. Using 3D relativistic hydrodynamical simulations, we explore the interaction of a relativistic jet, endowed with a luminosityL_j≈ 10^45.5erg s⁻¹and durationt_eng≈ 10 days compatible with those needed to power SLSNe, launched into the envelope of the exploding star. The jet successfully breaks through the expanding ejecta, and its shocked cocoon powers ultraviolet/optical emission lasting several days after the explosion and reaching a peak luminosity ≳10⁴⁴erg s⁻¹, corresponding to a sizable fraction ofL_j. This high radiative efficiency is the result of the modest adiabatic losses the cocoon experiences owing to the low optical depths of the enlarged ejecta at these late times, e.g., compared to the more compact stars in gamma-ray bursts. The luminosity and temperature of the cocoon emission match those of the “bumps” in SLSN light curves observed weeks prior to the optical maximum in many SLSNe. Confirmation of jet breakout signatures by future observations (e.g., days-long to weeks-long internal X-ray emission from the jet for on-axis observers, spectroscopy confirming large photosphere velocitiesv/c≳ 0.1, or detection of a radio afterglow) would offer strong evidence for central engines powering SLSNe. 

Abstract GRB 221009A is one of the brightest transients ever observed, with the highest peak gamma-ray flux for a gamma-ray burst (GRB). A Type Ic-BL supernova (SN), SN 2022xiw, was definitively detected in late-time JWST spectroscopy (t= 195 days, observer frame). However, photometric studies have found SN 2022xiw to be less luminous (10%−70%) than the canonical GRB-SN, SN 1998bw. We present late-time Hubble Space Telescope (HST)/WFC3 and JWST/NIRCam imaging of the afterglow and host galaxy of GRB 221009A att∼185, 277, and 345 days post-trigger. Our joint archival ground, HST, and JWST light-curve fits show strong support for a break in the light-curve decay slope att= 50 ± 10 days (observer frame) and a SN at <1.5× the optical/near-IR flux of SN 1998bw. This break is consistent with an interpretation as a jet break when requiring slow-cooling electrons in a wind medium with an electron energy spectral indexp> 2 andν_m<ν_c. Our light curves and joint HST/JWST spectral energy distribution (SED) also show evidence for the late-time emergence of a bluer component in addition to the fading afterglow and SN. We find consistency with the interpretations that this source is either a young, massive, low-metallicity star cluster or a scattered-light echo of the afterglow with a SED shape off_ν∝ν^2.0±1.0.